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The  Generation  of  Three-Dimensional  Body-Fitted 
Coordinate  Systems  For  Viscous  Flow  Problems 


by 

Z.  U.  A.  War si* 


Abstract 

An  analytical  model  for  the  generation  and  redistribution  of 
surface  coordinates,  which  is  to  be  used  along  with  the  full  3D  code, 
has  been  developed.  This  essentially  completes  the  development  of 
spatial  coordinate  generation  around  multibodies  and  particularly  around 
a  wing-body  combination  in  the  3D  space.  Numerical  results  for  some 
multibody  problems  have  been  obtained. 


Introduction 

The  problem  of  numerical  grid  generation  for  three-dimensional 
configurations  through  elliptic  PDE*s  has  been  pursued  under  this  grant 
for  the  last  several  years.  Various  reports  and  journal  publications 
have  explained  the  mathematical  model  in  detail.  This  model,  which  is 
based  on  the  numerical  solution  of  three  nonlinear  PDE’s  has  been  pro¬ 
gramed  on  Cray  1,  and  has  been  tested  for  single  and  two-body  configurations 
as  enclosed  in  a  single  outer  boundary. 

In  the  period  of  this  report  it  was  found  necessary  to  have  the 
capability  of  generating  those  surface  coordinate  systems  which  simplify 
the  procedure  of  3D  coordinates  generation  through  those  equations  which 
have  already  been  developed  in  the  previous  reporting  years.  The  main 
thrust  for  this  development  has  come  from  the  problem  of  wing-body 
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combination.  In  the  wing-body  configuration  we  must,  first  of  all,  generate 
a  suitable  surface  coordinate  system  which  provides  the  Dirichlet  boundary 
conditions  for  the  full  3D  code.  Therefore  a  coordinate  system  is  needed 
which  reflects  the  trace  of  the  body  and  wing  intersection  as  one  of  the 
coordinates,  and  the  coordinate  along  the  wing  as  the  outgoing  coordinate 
from  the  fuselage.  The  developed  analytical  model  which  has  been  described 
in  the  next  section  generates  those  surface  coordinates  on  the  fuselage 
in  which  the  trace  of  wing  and  body  intersection  is  one  of  the  surface 
curves.  Based  on  the  data  so  generated,  the  3D  equations  as  described 
in  References  4-8  have  been  solved.  Some  results  of  the  3D  coordinates 
for  various  multibodies  are  shown  in  Figures  1^5, 
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Equations  For  The  Transformed  Coordinates  In  A  Surface 


Let  (u,v)  be  a  coordinate  system  in  a  surface  (say  a  fuselage) 
on  which  the  coordinate  n  *  const.,  and  which  has  been  generated  by 
some  means,  e.g.,  from  Craidonfs  routine  (Ref.  9  ).  We  now  wish  to 
introduce  other  coordinates  £  *  £(u,v)  and  £  *  C(u,v)  in  the  same 
surface.  From  Warsi  (Ref.  4  ),  the  fundamental  equations  for  the 
generation  of  the  cartesian  coordinates  r  =  (x,  y,  z)  as  functions 
of  (C,  £),  are 


where 


Lr  +  Pr^  +  Qr^  =  nR  ,  n  «  const.,  (1) 


L  = 


n  = 


8liact  "  2g133C^  +  8333?C 
(X,  Y,  Z,),  surface  normal. 


Now 


811  =  U  *  U  ’  813  =  u  '  g33 =  u  * 

G2  =  gllg33  “  (g13* 


u 


r  u_  +  r  vr 

-u  c  -v  e 


» 


so  that  on  substituting  the  terms  r  #  r^,  r  in  (1)  and  collecting 
terms,  we  get  after  some  manipulations 
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(2) 


au«  '  2buU  +  CU£€  +  +  Q«c  "  4  V  * 

avcc  "  2bvcc  +  cv«  +  *e  +  Qv€  "  ^  V  *  (3) 


where 

J2  '  Vt  -  Vt  1 

and  a,b,c  are  the  coefficients  which  have  the  derivatives  ,  u^,  v^, 
v^t  also  x^,  x^,  etc.  and  A^u,  are  the  Belramians.  By  imposing  the 
appropriate  boundary  conditions  on  u  and  v  the  equations  (2)  and  (3) 
can  be  solved  as  a  2D  grid  generation  problem  for  the  new  coordinates 
£  and  This  routine  can  therefore  be  used  to  redistribute  the  surface 
coordinates  in  any  manner  from  any  given  surface  coordinates  u  and  v. 

The  preceding  equations  have  been  programed  as  a  subroutine  with 
the  general  3D  program  and  is  used  when  the  coordinates  for  a  wing- 
body  combination  are  to  be  generated.  The  trace  of  the  wing  with  the 
fuselage  is  considered  to  form  an  inner  closed  curve  and  the  outer  curve 
is  any  arbitrary  closed  curve  near  the  mid  section  of  the  fuselage  (Fig.  6), 
Equations  (2)  and  (3)  then  provide  the  curves  which  are  consistent  with 
the  3D  routine  for  the  generation  of  space  coordinates. 
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Pig.  1  Three-dimensional  coordinates  between  a  thin  ellipsoid  surrounded 
by  an  ellipsoid,  (a)  Front  section,  (b)  Meridlnal  section. 
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Three-d latent.  a«l  coordinates  between  a  thick  ellipsoid  surrounded 
by  an  ellipsoid,  (a)  Front  section,  (b)  Merldlnal  section. 
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two  Intersectlne  spheres. 


Fig.  4  Three-dimensional  coordinates  between  the  surface  of  Fig.  3  and 
an  outer  sphere,  (a)  Front  section  of  the  small  sphere, 

(b)  Front  section  at  the  Junction,  (c)  Front  section  of  the  large 
sphere,  (d)  Merldlnal  section. 


Fig.  6  Coordinate  redistribution  in  a  surface  obtained  by  solving  Eq 
and  (3).  (a)  Coordinates  ut  v  on  an  ellipsoid  through  Ref.  9 

(b)  Choice  of  the  inner  and  outer  curves  k  *  ^  and  £  *  ^2* 

(c)  Generated  coordinates  £  »  const.,  and  C  *  const. 


